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An ideal and reversible transfer technique for the quantum 
state between light and metastable collective states of matter 
is presented and analyzed in detail. The method is based on 
the control of photon propagation in coherently driven 3- level 
atomic media, in which the group velocity is adiabatically re- 
duced to zero. Form-stable coupled excitations of light and 
matter ( "dark-state polaritons" ) associated with the propaga- 
tion of quantum fields in Electromagnetically Induced Trans- 
parency are identified, their basic properties discussed and 
their application for quantum memories for light analyzed. 



I. INTRODUCTION 

Recent advances in quantum information science lead 
to many interesting new concepts such as quantum com- 
putation, quantum cryptography and teleportation . 
The practical implementation of quantum processing 
protocols requires coherent manipulation of a large num- 
ber of coupled quantum systems which is an extremely 
difficult task. One of the particular challenges for the 
implementation of these ideas involves physically trans- 
porting or communicating quantum states between dif- 
ferent nodes of quantum networks 0] . Quantum optical 
systems appear to be very attractive for the realization 
of such networks. On one hand photons are ideal car- 
riers of quantum information: they are fast, robust and 
readily available. On the other hand atoms represent re- 
liable and long-lived storage and processing units. There- 
fore the challenge is to develop a technique for coherent 
transfer of quantum information carried by light to atoms 
and vise versa. In other words it is necessary to have a 
quantum memory that is capable of storing and releasing 
quantum states on the level of individual qubits and on 
demand. Such a device needs to be entirely coherent and 
in order to achieve a unidirectional transfer (from field to 
atoms or vise versa) an explicit time dependent control 
mechanism is required. 

Classical optical data storage in the time domain, 
based on the phenomenon of spin- j^] and photon echo Q , 
has a long history. After the first proposals of stimulated 
two-level photon echo JjJ and demonstrations of light- 
pulse storage in these systems H many important devel- 
opments have taken place in this field. Particularly inter- 
esting are techniques based on Raman photon echos ft as 



they combine the long livetime of ground-state hyperfine 
or Zeeman coherences for storage with data transfer by 
light at optical frequencies |l(|. While these techniques 
are very powerful for high-capacity storage of classical 
optical data, they cannot be used for quantum memory 
purposes. The techniques employ direct or dressed-state 
optical pumping and thus contain dissipative elements 
or have other limitations in the transfer process between 
light and matter. As a consequence they do not operate 
on the level of individual photons and cannot be applied 
for quantum information processing. 

The conceptually simplest approach to a quantum 
memory for light is to "store" the state of a single pho- 
ton in an individual atom. This approach involves a co- 
herent absorption and emission of single photons by sin- 
gle atoms. However, the single-atom absorption cross- 
section is very small, which makes such a process very 
inefficient. A very elegant solution to this problem is 
provided by cavity QED . Placing an atom in a high- 
er resonator effectively enhances its cross-section by the 
number of photon round-trips during the ring-down time 
and thus makes an effective transfer possible. Raman 
adiabatic passage techniques p2| with time-dependent 
external control fields can be used to implement a di- 
rected but reversible transfer of the quantum state of a 
photon to the atom (i.e. coherent absorption). However, 
despite the enormous experimental progress in this field 
p3| , it is technically very challenging to achieve the nec- 
essary strong-coupling regime. Furthermore the single- 
atom system is by construction highly susceptible to the 
loss of atoms and the speed of operations is limited by 
the large Q-i actor. 

On the other hand a photon can be absorbed with 
unit probability in an optically thick ensemble of atoms. 
Normally such absorption is accompanied by dissipative 
processes which result in decoherence and thus deterio- 
rate the quantum state. Nevertheless it has been shown 
that such absorption of light leads to a partial mapping 
of its quantum properties to atomic ensembles Jl4] , p^| . 
As a consequence of dissipation these methods do not al- 
low to reversibly store the quantum state on the level of 
individual photon wave-packets (single qubits). Rather 
a stationary source of identical copies is required (e.g. a 
stationary source of squeezed vacuum, which can be con- 
sidered as a train of identical wave-packets in a squeezed 
vacuum state) to partially map quantum statistics from 
light to matter. 

Recently we have proposed a method that combines 
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the enhancement of the absorption cross section in many- 
atom systems with dissipation-free adiabatic passage 
techniques ]T^-[l8t. It is based on an adiabatic transfer 
of the quantum state of photons to collective atomic ex- 
citations using electromagnetically induced transparency 
(EIT) in 3- level atoms [jl9). Since the technique alle- 
viates most of the stringent requirements of single-atom 
cavity-QED, it could become the basis for a fast and reli- 
able quantum network. Recent experiments [ pO|pl| ] have 
already demonstrated one of the basic principle of this 
technique - the dynamic group velocity reduction and 
adiabatic following in the so-called "dark-state" polari- 
tons. The aim of the present and subsequent papers is 
to analyze the physics of the reversible storage technique 
in detail and to discuss its potentials and limitations. 

Electromagnetically induced transparency can be used 
to make a resonant, opaque medium transparent by 
means of quantum interference. Associated with the 
transparency is a large linear dispersion, which has been 
demonstrated to lead to a substantial reduction of the 
group velocity of light Q . Since the group velocity re- 
duction is a linear process, the quantum state of a slowed 
light pulse can be preserved. Therefore a non-absorbing 
medium with a slow group velocity is in fact a tempo- 
rary "storage" device. However, such a system has only 
limited "storage" capabilities. In particular the achiev- 
able ratio of storage time to pulse length is limited by 
the square root of the medium opacity p3[ and can prac- 
tically attain only values on the order of 10 to 100. This 
limitation originates from the fact that a small group 
velocity is associated with a narrow spectral acceptance 
window of EIT p4| and hence larger delay times require 
larger initial pulse length. 

The physics of the state-preserving slow light propa- 
gation in EIT is associated with the existence of quasi- 
particles, which we call dark-state polaritons (DSP). A 
dark-state polariton is a mixture of electromagnetic and 
collective atomic excitations of spin transitions (spin- 
wave). The mixing angle between the two components 
determines the propagation velocity and is governed by 
the atomic density and the strength of an external con- 
trol field. The key idea of the present approach is the 
dynamic rotation of the mixing angle which leads to an 
adiabatic passage from a pure photon-like to a pure spin- 
wave polariton thereby decelerating the initial photon 
wavepacket to a full stop. In this process the quantum 
state of the optical field is completely transferred to the 
atoms. During the adiabatic slowing the spectrum of the 
pulse becomes narrower in proportion to the group ve- 
locity, which essentially eliminates the limitations on ini- 
tial spectral width or pulse length and very large ratios 
of storage time to initial pulse length can be achieved. 
Reversing the rotation at a later time regenerates the 
photon wave-packet. Hence the extension of EIT to a 
dynamic group-velocity reduction via adiabatic follow- 
ing in polaritons can be used as the basis of an effective 



quantum memory. Before proceeding we note some ear- 
lier work on the subject. The polariton picture of Raman 
adiabatic passage has first been introduced in Ref. p5| . 
Furthermore Grobe and coworkers [ p6[ pointed out that 
the spatial profile of an atomic Raman coherence can be 
mirrored into the electromagnetic field by coherent scat- 
tering, whereas time-varying fields can be used to create 
spatially non-homogeneous matter excitations. 

In the present paper we will present a quantum picture 
of slow light propagation in EIT in terms of dark-state 
polaritons. We will analyze the properties of the po- 
laritons and discuss their application to reversible, fast 
and high-fidelity quantum memories. Limitations and re- 
strictions of the transfer process from non-adiabatic pro- 
cesses will be discussed as well as effects from the medium 
boundary and atomic motion. Other important aspects 
of the collective quantum memory such as its decoherence 
properties will be the subject of subsequent publications. 



II. QUANTUM MEMORY FOR A SINGLE-MODE 
FIELD 

The essential aspects of the quantum state mapping 
technique can be most easily understood for the case of 
a single mode of the radiation field as realized e.g. in 
a single mode optical cavity. In what follows we will 
address this case first in order to motivate the following 
discussion on propagating photon wave packets. 

Consider a collection of N 3-level atoms with two meta- 
stable lower states as shown in Fig. |l| interacting with 
two single-mode optical fields. The transition \a) — ► \b) 
of each of these atoms is coupled to a quantized radi- 
ation mode. Moreover the transitions from |a) — * |c) 
are resonantly driven by a classical control field of Rabi- 
frequency Q. The dynamics of this system is described 
by the interaction Hamiltonian: 



N 

H = hg^2aa l a 

i=l 



N 



h.c. 



(1) 



Here cr* = \n)a(v\ is the flip operator of the ith atom 
between states \fi) and \u). g is the coupling constant 
between the atoms and the quantized field mode (vacuum 
Rabi-frequency) which for simplicity is assumed to be 
equal for all atoms. 

When all atoms are prepared initially in level \b) the 
only states coupled by the interaction are the totally sym- 
metric Dicke-like states [E7I 
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In particular, if the field is initially in a state with at most 
one photon, the relevant eigenstates of the bare system 
are: the total ground state |b, 0), which is not affected by 
the interaction at all, the ground state with one photon 
in the field |b, 1), as well as the singly excited states |a, 0) 
and |c, 0). For the case of two excitations, the interaction 
involves 3 more states etc. The coupling of the singly and 
doubly excited systems are shown in Fig. [l] b. 



a) 



b) 




|b,0> 




ac,0> 



FIG. 1. a) 3-level atoms coupled to single quantized mode 
and classical control field of (real) Rabi-frequency Q(t). b) 
coupling of relevant bare eigenstates for at most one photon 

The interaction has families of dark-states, i.e. states 
with zero adiabatic eigenvalue [ pdp^ , p^ ]. The simplest 
one is 

|Z>,1) =cos0(t)|b,l) -sin0(i)|c,O), (6) 



tan 0(f) 



and in general one has 
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The dark states do not contain the excited state and are 
thus immune to spontaneous emission. It should also be 
noted that although the dark states \D, n) are degenerate 
they belong to exactly decoupled subsystems as long as 
spontaneous emission is disregarded. This means there 
is no transition between them even if non-adiabatic cor- 
rections are taken into account. The existence of collec- 
tive dark states provides a very elegant way to transfer 
the quantum state of the single-mode field to collective 
atomic excitations. Adiabatically rotating the mixing an- 
gle from to tt/2 leads to a complete and reversible 
transfer of the photonic state to a collective atomic state 
if the total number of excitations n is less than the num- 
ber of atoms. This can be seen very easily from the ex- 
pression for the dark states, eq. (|8]) : If 8 : — > tt/2 one 
has for all n < N 



Thus if the initial quantum state of the single- mode light 
field is in any mixed state described by a density matrix 
Pf = J2nmP nm l n )( m l> the transfer process generates a 
quantum state of collective excitations according to 
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It should be noted that the quantum-state transfer does 
not necessarily constitute a transfer of energy from the 
quantum field to the atomic ensemble. Since in the Ra- 
man process the coherent "absorption" of a photon from 
the quantized mode is followed by a stimulated emission 
into the classical control field, most of the energy is ac- 
tually deposited in the latter field. 

The transfer of quantum states between light and mat- 
ter due to adiabatic following in collective dark states is 
the key point of the present work. Before proceeding we 
also note that the transfer rate is proportional to the to- 
tal number of atoms N, which is a signature of collective 
coupling. This makes the proposed method potentially 
fast and robust. 



III. QUANTUM DESCRIPTION OF 
SLOW-LIGHT PROPAGATION 



We now discuss a generalization of the mapping tech- 
nique to propagating fields. The adiabatic transfer of 
the quantum state from the radiation mode to collec- 
tive atomic excitations discussed in the previous section 
is strongly related to intracavity electromagnetically in- 
duced transparency (EIT) p8fl . In order to generalize the 
technique to multi-mode fields it is useful to discuss first 
the propagation of light in 3-level media under conditions 
of EIT. 



A. model 

Consider the quasi 1-dimensional problem shown in 
Fig. |^. A quantized electromagnetic field with positive 
frequency part of the electric component JS'+) couples 
resonantly the transition between the ground state \b) 
and the excited state |<x). v = uj a b is the carrier fre- 
quency of the optical field. The upper level \a) is fur- 
thermore coupled to the stable state |c) via a coherent 
control field with Rabi-frequency f2. 



3 




(CO - COab ) / Y 

FIG. 2. top: 3-level A-type medium resonantly coupled to 
a classical field with Rabi-frequency Q(t) and quantum field 
E(z,t). bottom: Typical susceptibility spectrum for probe 
field E as function of normalized detuning from resonance for 
resonant drive field. Real part x' describes refractive index 
contribution and imaginary part x" absorption. 

The interaction Hamiltonian reads 

v = -pJ2( a i b ^ +) ^) + h - a ) 

-h^facWzj^e^-^+h.a) (11) 

3 

where zj denotes the position of the jth atom, p denotes 
the dipole matrix element between the states \a) and |6), 
and 
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defines the atomic flip operators, k d — kd ■ e z — cos$ 
is the projection of the wavevector of the control field to 
the propagation axis of the quantum field. For the sake 
of simplicity we here assume that the carrier frequencies 
v and Ud of the quantum and control fields coincide with 
the atomic resonances tu ab and u) ac respectively. Motional 
effects and the associated Dopplcr-shifts will be discussed 
later. We introduce slowly-varying variables according to 
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Here V is some quantization volume, which for simplicity 
was chosen to be equal to the interaction volume. 

If the (slowly-varying) quantum amplitude does not 
change in a length interval Az which contains N z 3> 1 
atoms we can introduce continuum atomic variables 



iVz z jf =N, 



and make the replacement X^jLi — > T I wnere N 
is the number of atoms, and L the length of the inter- 
action volume in propagation direction of the quantized 
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field. This yields the continuous form of the interaction 
Hamiltonian 

V = - J^(hgNa ab (z,t)£(z,t) 

+Ml(z, t)c lAkz N(z)a ac (z, t) + h.a) . (16) 



Here g = P^J 2 he v ' 1S tnc a t° m - nc ld coupling constant 

and Afc = fcj - k d = ^f- (cos $ — 1). 

The evolution of the Heisenberg operator correspond- 
ing to the quantum field can be described in slowly vary- 
ing amplitude approximation by the propagation equa- 
tion 

d d \ 

g- t +c—j£(z,t) = igNa ba (z,t). (17) 

The atomic evolution is governed by a set of Heisenberg- 
Langevin equations 

-i(n* e - iAkz a ca -h.aj+F a , (18) 
&bb = ja aa + ig(s f a ba - h.a^j + F b , (19) 
°cc = i° aa + i (n*e- tAkz a ca - /i.a) + F c , (20) 

Cfca = —"fbaVha + ig£ (&bb — &aaj 

+ifle lAkz a bc + F ba , (21) 



A.. . — 't~\ iAhz ( ' — ■ — \ 

0"ca = -Ica^ca + " 2C |cr cc -(T aa l 

a bc = in*e~ zAkz a ba - ig£<r ac , 



(22) 
(23) 



7 a = 7 + 7' and 7,7' denote longitudinal and 7 M „ 
transversal decay rates. and F^ u are e)-correlated 
Langevin noise operators, whose explicit form is not of 
interest here. 

It should be noted that we have disregarded dissipative 
population exchange processes due to e.g. spin-Hip col- 
lisions and dephasing of the lower-level transition. This 
is justified since we assume that the interaction time is 
sufficiently short compared to the characteristic times of 
these processes. Both, dissipative population exchange 
as well as dephasing will be discussed in detail later on. 



B. low-intensity approximation 

In order to solve the propagation problem, we now as- 
sume that the Rabi-frequency of the quantum field is 
much smaller than Q and that the number density of 
photons in the input pulse is much less than the number 
density of atoms. In such a case the atomic equations 
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can be treated perturbatively in £. In zeroth order only 
a bb = 1 is different from zero and in first order one finds 



iAk 



. z d_ 
dt 



0~ be 



(24) 



With this the interaction of the probe pulse with the 
medium can be described by the amplitude of the probe 
electric field £ and the collective ground-state spin vari- 
able a bc : 



^ + ^ ] E(z,t) = §e^ bc 
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dz 



(25) 
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C. adiabatic limit 

The propagation equations simplify considerably if we 
assume a sufficiently slow change of fi, i.e. adiabatic 
conditions [^9 23 3(J. Normalizing the time to a char- 
acteristic scale T via t = t/T and expanding the r.h.s. 
of (|2^) in powers of 1/T we find in lowest non- vanishing 
order 



a bc {z,t) = -5^ e 



(27) 



We note that also the noise operator F ba gives no con- 
tribution in the adiabatic limit, since (F x (t)F y (t')) ~ 
6{t - t') = 5(t- t') jT. Thus in the perturbative and adi- 
abatic limit the propagation of the quantum light pulse 
is governed by 
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0i + C 8- z ^ M -- 



g 2 N d £( z ,t) 
n* dt n 



(28) 



D. slow light and delay-time limitations 

If Q(z,t) — f2(z) is constant in time, the term on the 
r.h.s. of the propagation equation ( pq ) simply leads to a 
modification of the group velocity of the quantum field 
according to 



n g (z) 
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(29) 
(30) 



with p being the atom density and A the resonant wave- 
length of the a — > b transition. The solution of the wave 
equation is in this case 



£{z,t) =£[0,t 



r dz' 

l0 Vg{z') 



(31) 



where £(0,t') denotes the field entering the interaction 
region at z = 0. This solution describes a propagation 
with a spatially varying velocity v g . It is apparent that 
the temporal profile of the pulse is unaffected by the slow- 
down. As a consequence the integrated electromagnetic 
energy flux through a plane perpendicular to the prop- 
agation is the same at any position. Furthermore the 
spectrum of the pulse remains unchanged 

S(*,w)= f Are- luJT (&(z,t)£(z, t - t)\ 

= S(0,w). (32) 
In particular the spectral width stays constant 



Auj p (z) = Acjp(0). 



(33) 



On the other hand a spatial change of the group velocity, 
either in a stepwise fashion as e.g. at the entrance of the 
medium or in a continuous way e.g. due to a spatially 
decreasing control field, leads to a compression of the 
spatial pulse profile. In particular, if the group velocity 
is statically reduced to a value v g , the spatial pulse length 
Al is modified according to 
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(34) 



as compared with a free-space value AIq. 

It is instructive to discuss the limitations to the achiev- 
able delay (and therefore "storage") time Td in an EIT 
medium with a very small but finite group velocity. A 
first and obvious limitation is set by the finite lifetime of 
the dark state, which had been neglected in the previ- 
ous subsection. If j bc denotes the dephasing rate of the 
b — c transition it is required that < r y^ c 1 ■ A much 
stronger limitation arises however from the violation of 
the adiabatic approximation. 

The EIT medium behaves like a non-absorbing, lin- 
ear dispersive medium only within a certain frequency 
window around the two-photon resonance (see Fig.^ and 
[B4J ). The adiabatic approximation made in the last 
sub-section essentially assumes that everything happens 
within this frequency window. If the pulse becomes too 
short, or its spectrum too broad relative to the trans- 
parency width, absorption and higher-order dispersion 
need to be taken into account. 

The transparency window is defined by the intensity 
transmission of the medium. In order to determine its 
width we consider the susceptibility of an ideal, homoge- 
neous EIT medium with a resonant drive field. Here one 
has 
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(35) 



where S = v — u> ac is the detuning of the probe field. If 
we assume a homogeneous drive field, we find 



T(5, z) = cxp| — fczlm[x]| 
« exp{-5 2 /Au£}, 



with 



Au>tr = 
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I being the propagation length in the medium and a = 
■J^zp\ 3 kl the opacity in the absence of EIT. One recog- 
nizes that the transparency width decreases with increas- 
ing group index. It is instructive to express the trans- 
parency width in terms of the pulse delay time = n g l/c 
for the medium. This yields 
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(38) 



Hence large delay times imply a narrow transparency 
window, which in turn requires a long pulse time. When 
the group velocity gets too small such that the trans- 
parency window becomes smaller than the spectral width 
of the pulse, the adiabatic condition is violated, and the 
pulse is absorbed. Hence there is an upper bound for 
the ratio of achievable delay (storage) time to the initial 
pulse length of a photon wavepacket 



— <\[a and r d < jZ 



(39) 



The ratio r^/Vp is the figure of merit for a memory device. 
The larger this ratio the better suited is the system for 
a temporal storage. 

There is another quantity which is important for the 
storage capacity, namely the ratio of medium length L to 
the length L p of an individual pulse inside the medium. 
Following similar arguments as above one finds that this 
quantity is also limited by the square root of the opacity 

^ < yfa, (40) 

Lip 

where n g 3> 1 was assumed. 

In practice, the achievable opacity a of atomic vapor 
systems is limited to values below 10 4 resulting in upper 
bounds for the ratio of time delay to pulse length of the 
order of 100. Thus dense EIT media with ultra-small 
group velocity are only of limited use as a temporary 
storage device. The propagation velocity cannot be made 



zero and non-adiabatic corrections limit the achievable 
ratio of storage time to the time length of an individual 
qubit. 

It should be mentioned that the narrowing of the EIT 
transparency window is also a consideration for effects in- 
volving freezing light in a moving media [ |32| and so-called 
optical black holes based on EIT |33| . In the following we 
will show, that EIT can nevertheless be used for mem- 
ory purposes in a very effective way when combined with 
adiabatic passage techniques. 



IV. DARK-STATE POLARITONS 

In the previous section we have discussed the propa- 
gation of a quantum field in an EIT medium under oth- 
erwise stationary conditions, i.e. with a constant or only 
spatially varying control field. Since under these condi- 
tions, the hamiltonian of the system is time-independent, 
a coherent process that allows for an uni-directional 
transfer of the quantum state of a photon wavepacket 
to the atomic ensemble was not possible. We will show 
now that this limitation can be overcome easily by allow- 
ing for a time- dependent control field. This provides an 
elegant tool to control the propagation of a quantum light 
pulse. For a spatially homogeneous but time-dependent 
control field, f2 = £l(t), the propagation problem can 
be solved in a very instructive way in a quasi-particlc 
picture. In the following we will introduce these quasi- 
particles, called dark-state polaritons 0,^1 and discuss 
their properties, applications and limitations. 



A. definition of dark- and bright-state polaritons 

Let us consider the case of a time-dependent, spatially 
homogeneous and real control field fl — £l(t) = il(t)* . 
We introduce a rotation in the space of physically rele- 
vant variables - the electric field £ and the atomic spin 
coherence <7b c - defining two new quantum fields ^(z,t) 
and $(z, t) 



* = cos 9(t) £(z, t) - sin 6(t) VNa bc (z, t) e lAkz , (41) 
$ = sin0(i) S(z, t) + cos 9(t) ViV5 6c (z, t) e lAkz , 

with the mixing angle 9(t) given by 



(42) 



tan 2 9(t) = 



g 2 N 

n 2 {t) 



(43) 



i& and <t are superpositions of electromagnetic (£) and 
collective atomic components (\/N(Tbc), whose admixture 
can be controlled through 9(t) by changing the strength 
of the external driving field. 

Introducing a plain-wave decomposition ^{z,t) = 
T, k ^k(t)e ikz and $(z,t) = Y. k ®k\t)e lkz respectively, 
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one finds that the mode operators obey the commuta- 
tion relations 



B. adiabatic limit 



cos 



= &k,k 

= 5k,k 
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(44) 
(45) 
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In the linear limit considered here, where the number 
density of photons is much smaller than the density of 



atoms, <x 



1 <7 J 



0. Thus the new fields possess 



bosonic commutation relations 



0, 



(47) 
(48) 



and we can associate with them bosonic quasi-particles 
(polaritons). Furthermore one immediately verifies that 
all number states created by 
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(49) 



where |0) denotes the field vacuum, are dark-states 
PUP^f . The states \rik) do not contain the excited atomic 
state and are thus immune to spontaneous emission. 
Moreover, they are eigenstates of the interaction Hamil- 
tonian with eigenvalue zero, 



V\n k )=0. 



(50) 



For these reasons we call the quasi-particles \& "dark- 
state polaritons" . Similarly one finds that the elementary 
excitations of $ correspond to the bright-states in 3-level 
systems. Consequently these quasi-particles are called 
"bright-state polaritons" . 

One can transform the equations of motion for the elec- 
tric field and the atomic variables into the new field vari- 
ables. In the low-intensity approximation one finds 
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ba : 



where one has to keep in mind that the mixing angle 9 is 
a function of time. 



Introducing the adiabaticity pa- 

rameter e = (q^/NT^ with T being a characteristic 

time, one can expand the equations of motion in powers 
of e. In lowest order i.e. in the adiabatic limit one finds 



$ w 0. 

Consequently 

S(z,t)= cos6(t)4>(z,t), 



(53) 



(54) 

'Nabc = -sm9(t)4>(z,t)e- tAkz . (55) 
Furthermore \& obeys the very simple equation of motion 

— +ccos 2 9(t) — *(2,t)=0. (56) 



C. "stopping" and re-accelerating photon 
wavepackets 

Eq. ( [56"1) describes a shape- and quantum-state preserv- 
ing propagation with instantaneous velocity v = v g (t) = 
ccos 2 9(t): 



drcos 2 6»(r),0 



(57) 



For 9^0, i.e. for a strong external drive field f2 2 3> 
g 2 N, the polariton has purely photonic character ^ = 
£ and the propagation velocity is that of the vacuum 
speed of light. In the opposite limit of a weak drive field 
SI 2 <C g 2 N such that 9 — ► ir/2, the polariton becomes 
spin- wave like $ = —\ZNab c e lAkz and its propagation 
velocity approaches zero. Thus the following mapping 
can be realized 



£{z) <=> a bc (z')e 



■Akz 



(58) 



with z 1 = z + z a = z + J °°drccos 2 9(t). This is the 
essence of the transfer technique of quantum states from 
photon wave-packets propagating at the speed of light 
to stationary atomic excitations (stationary spin waves). 
Adiabatically rotating the mixing angle from 9 = to 9 = 
ir/2 decelerates the polariton to a full stop, changing its 
character from purely electromagnetic to purely atomic. 
Due to the linearity of eq.(|56|) and the conservation of 
the spatial shape, the quantum state of the polariton is 
not changed during this process. 

Likewise the polariton can be re-accelerated to the vac- 
uum speed of light; in this process the stored quantum 
state is transferred back to the field. This is illustrated 
in Fig.^J, where we have shown the coherent amplitude 
of a dark-state polariton which results from an initial 
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light pulse as well as the corresponding field and matter 
components. 



0.5 



~A e 
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50 { 100 
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FIG. 3. Propagation of a dark-state polariton with 
envelope exp{ — (jz/10) 2 }. The mixing angle is rotated 
from to 7r/2 and back according to cot 9(t) — 
100(1 - 0.5 tanh[0.1(t- 15)] +0.5tanh[0.1(i- 125)]) as shown 
in (a). The coherent amplitude of the polariton ^ = is 
plotted in (b) and the electric field E = (E) and matter com- 
ponents \a c t\ = \{&cb)\ in (c) and (d) respectively. Axes are 
in arbitrary units with c = 1. 



D. simultaneous narrowing of transparency window 
and pulse spectral width 

As was discussed above, the transparency window of 
the EIT medium, i.e. the range of frequencies for which 
absorption is negligible, decreases with the group veloc- 
ity. Since the bandwidth of the propagating pulse should 
always be contained within this range of frequencies to 
avoid absorption, the question arises whether this pre- 
vents the stopping of the polariton. 

To answer this question we first note that during the 
process of adiabatic slowing the spatial profile, and in 
particular the length of the wavepacket (A/) remains un- 
affected, as long as the group velocity is only a function 
oft. I.e. 



a; = a/ . 



(59) 



At the same time, the amplitude of the electric field gets 
reduced and its temporal profile stretched due to the 
reduction of the group velocity. The opposite happens 
when the group velocity is increased. One finds from 



£(z,t) 



cos 6(t) 
cos (9(0) 



£ [ z — c dr cos 2 



.0 



(60) 



As a consequence the spectrum of the probe field changes 
during propagation. Assuming that cos0 changes only 
slow compared to the field amplitude one finds 

cos 2 6{t) 



S{z,uj) 



cos 2 0(0) 



S 0, 



UJ 



cos 2 6{t) 



(61) 



In particular the spectral width narrows (broadens) ac- 
cording to 



Awp(t) w Awp(0) 



cos 2 d(t) 
cos 2 0(0) ' 



(62) 



When the group velocity of the polariton is reduced in 
time due to reduction of the control-field amplitude the 
EIT transparency window shrinks according to 



Auj t r(t) 



cot 2 0(f) 
cot 2 0(0) 



Aw tr (0). 



(63) 



However the spectral width of the wavepacket shrinks as 
well and the ratio remains finite: 



Awp(i) _ sin 2 0(7) Awp(0) 



: 0(O) Aw tr (0)' 



(64) 



As will be discussed later on, for practically relevant 
cases, sin 2 9(t)/ sin 2 0(0) is always close to unity. Thus 
absorption can be prevented in the dynamic light- 
trapping method as long as the input pulse spectrum 
lies within the initial transparency window: 



Awp(0) < Aw tr (0). 



(65) 



As we have already noted earlier this condition can easily 
be fulfilled if an optically dense medium is used. The 
simultaneous reduction in transparency bandwidth and 
pulse bandwidth is illustrated in Fig.^J. 



7-77 CO ) 



S((0) 




CO / co 

FIG. 4. Simultaneous narrowing of transmission spectrum 
(top) and pulse spectrum (bottom) for time dependent varia- 
tion of group velocity v/c in units of u>o = g 2 N/j. Parameters 
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E. boundary behavior and initial pulse compression 

In the above discussion we have analyzed the propa- 
gation of the probe pulse inside the medium. We now 
turn to the behavior at the medium boundaries. There 
are two issues of interest: (i) reflection from the medium 
surface and (ii) the effects of a possible steplike change 
of the group velocity when passing from vacuum or air 
to the medium. 

Under ideal conditions the refractive index of the EIT 
medium is exactly unity on resonance and hence there 
is no reflection for the resonant component of the in- 
put pulse. If the spectrum of the input pulse is further- 
more sufficiently narrow as compared to the initial trans- 
parency window of EIT, the refractive index is very close 
to unity over the entire relevant bandwidth and no field 
component gets reflected. To quantify this we note, that 
the index of refraction near resonance of an idealized, 
resonantly driven 3-level medium can be written as 



n{oj) 




(66) 



The reflection coefficient for normal incidence from vac- 
uum to the medium given by 



is plotted in fig||. 



1 — n{uj) 



1 + n{uj) 



(67) 



^ 0.4 

0.2 
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FIG. 5. Reflection coefficient for normal incidence as func- 
tion of probe detuning from resonance. 

Near resonance R(u>) can be approximated by 



( v° 



Auj 



(68) 



where Au> = u — uj a b- Since typically Aw p (0) <C uj a bVg/c 
reflection can be neglected. 

When the group velocity inside the medium is smaller 
than c at the time when the probe pulse enters the 
medium, the pulse gets spatially compressed according 
to eq.(M). This substantially reduces the requirements 



on the size of the medium and enhances the storage ca- 
pacity. 

Continuity of the electric field amplitude £ at the en- 
trance of the medium (z = 0) implies a jump in the 
polariton amplitude 



$(+o,t) = /4*HM) 

V v% 



(69) 



Hence the density of polaritons increases by a factor c/v g 
when entering the medium, while there total number is 
conserved since the pulse length is changed by the inverse 
factor. 

Boundary effects can also be used for a controlled de- 
formation of a stored light pulse. When the polariton 
reaches the back end of the medium before cos 0(f) ap- 
proaches unity, the pulse will be spatially modulated. 



V. LIMITATIONS OF QUANTUM 
STATE-TRANSFER 

In the preceeding sections we have discussed the con- 
trol of propagation of the photon wavepackets under the 
assumptions of an adiabatic evolution and small probe- 
field amplitudes. We furthermore neglected ground- 
state decoherence and motional (Doppler) effects. In the 
present section we will discuss the validity of these ap- 
proximations in detail. 



A. non-adiabatic corrections 

Let us first discuss the effect of non-adiabatic transi- 
tions i.e. let us take into account terms in first order of e 
in eqs.(|5l|) and (p2|) . For the following discussion it is suf- 
ficient to consider the semiclassical limit, where the op- 
erator character of the variables as well as the Langevin 
noise Fb a can be disregarded. Up to e 1 one obtains 



$ ps sin 2 0(i)^-0* 
9 2 N 



sin 3 0(t) 7 
cos0(t) ~c?N 



(70) 



On the same level of approximations we can replace \E' 
on the right hand side by — ccos 2 Od^S/dz. This gives 



$ fa sin 2 0(t) 



- a "> 7 ^-sin 3 0(t)cos0(i)^-c— - (71) 



9 2 N 



g 2 N dz 



Substituting this into eq. (^fj) yields the lowest-order non- 
adiabatic corrections to the propagation equation of the 
dark-state polariton: 



— +ccos 2 (t — 
at oz 



r> 2 B 3 



(72) 
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where 



A(t) 

B(t) 
C(t) 
D(t) 



1 d \ 
1+ 2dl) 



¥ sin 2 i 



9 2 N 



sin 6 d 2 



sin J 8, 



3g 2 N dt 2 

1 d \ sin 4 cos 2 8 



g 2 N 



sin 4 8 cos 4 i 



(73) 

(74) 
(75) 
(76) 



A(t) describes homogeneous losses due to non-adiabatic 
transitions followed by spontaneous emission. B(t) gives 
rise to a correction of the polariton propagation veloc- 
ity, C(t) results in a pulse spreading by dissipation of 
high spatial- frequency components and D(t) leads to a 
deformation of the polariton. 

Since all coefficients depend only on time, the propa- 
gation equation can be solved by a Fourier transform in 
space *(z,t) = / dk^(k,t)e~ ikz . This yields 



®(k,t) = ^(k,0) expjifc J dt' v gl .(t')+cB(t') 
x expi-ik 3 c 3 Jdt'D(t')^ x 



(77) 



x exp 



dt' 



A(t') + k 2 c 2 C(t') 



The last term contains all losses due to non-adiabatic 
corrections. In order to neglect dissipation, the integral 
in the exponent of this term needs to be small compared 
to unity. Taking into account 8(0) = 8(oo) = this 
results in the two conditions 



7 fc 2 c 2 



° sin* 9(t) cos* 6(t) 

g 2 N « L ' 



and 



7 / dt 
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9 2 sin 2 ( 
9 2 N 



dt- 



o g 2 N + W(t) 



< 1. 



(78) 



(79) 



The first condition can be brought into a transparent 
form when replacing sin 4 9 by unity, which results in 



2„2 



7 Arc 



dt 



cos 2 8(t) jk 2 cz 



g 2 N 



g 2 N 



< l 



(80) 



for all relevant k. Here z = c J °° dt cos 2 8(t) is the prop- 
agation depth of the polariton inside the medium, and we 
have assumed cos 8(0) = 1 and cos#(oo) = for simplic- 
ity. Noting that the range of relevant spatial Fourier fre- 
quencies is determined by the inverse of the initial pulse 
length L p , i.e. setting k 



L p , one finds 



9 2 N t2 
z < Lt 

7c F 



or z <C VaLp, 



(81) 



where a = g 2 Nz/-fc is again the opacity of the medium 
without EIT. To illustrate this limitation we have shown 
in Figjf] the slow-down and successive acceleration of a 
light pulse obtained from a numerical solution of the full 
1-d propagation problem as well as the analytical ap- 
proximation following from eq. (|77j) . Here ^-^fj-^- = 2 
(z max = 250). One clearly recognizes a decay of the elm. 
excitation E long before the deceleration and the asso- 
ciated transfer to the atom systems sets in. The polari- 
ton energy 1$, = J dz\^(z,t)\ 2 decreases approximately 
exponentially with the propagation distance Az of the 
pulse center. In the example 1^ decays after Az = 250 
approximately by a factor e -0 322 ps 0.724. 




250 300 



FIG. 6. Deceleration (top) and re-aceleration (bottom) 
of a Gaussian light pulse exp[-(t - 500) 2 /(200) 2 ] with 
cot0(t) = 0.363(l-(2/7r)arccot[5(l-0.5tanh[0.005(t-2000)] 
+0.5tanh[0.005 * (t - 3200)])]). The parameters are 
g 2 N/'y 2 = 2.5, a = 625, c = 1. Pulses are shown for t =600, 
800, 1000, 1200, 1400, 1600, 1800, 2000, 2200, and 2400 in the 
top figure and for t =2800, 3000, 3200, 3400, 3800, 4000 in 
the bottom figure. One recognizes good agreement between 
numerical solution of the propagation equations (solid line) 
and analytic approximation (dashed line). 

Since the spatial pulse length in the medium is related 
to the initial group velocity v® (at time t = before the 
trapping) and to the initial pulse bandwidth as L p ~ 
Vg/Au! p (0) and since Aw tr (0) = \fajTd ~ yfav®/z, for 



c, equation (|81|) is simply another way to state: 



Au p < Aw tr (0). 



(82) 



Hence the present analysis confirms the qualitative con- 
ditions for adiabatic following: the initial probe spectrum 
should be contained within the original transparency 
window. Once again, this can be satisfied if the medium 
is optically dense, a ^>> 1. 

There is also a second condition ( 179 ) which is well- 
known from adiabatic passage requirement p9,p4| and 
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sets a limit to the rotation velocity 9 of the mixing angle 
and hence to the deceleration/acceleration of the polari- 
ton. Introducing a characteristic time scale of the accel- 
eration/deceleration period T we obtain 

T>^^, (83) 
c c 

where l a bs — crf/g 2 N is the absorption length in the ab- 
sence of EIT. We note that for realistic experimental pa- 
rameters the quantity on the right hand side is extremely 
small on all relevant time scales. Hence, in practice the 
rate of change of the mixing angle does not significantly 
affect the adiabatic dynamics. 

B. weak field approximation 

The analysis of sec. IV also involves a perturbation ex- 
pansion, valid when the control field is much stronger 
than the probe field. It is easy to see that this ex- 
pansion is justified even when the control-field Rabi- 
frequency is reduced to zero. Making use of (g7j) one 
finds: g 2 (E + E) / \fl\ 2 — (a cb a bc ) — (a cc ). I.e., the ratio 
of the average intensities of quantum and control field is 
proportional to the probability to find an atom in state c. 
If the initial number-density of photons in the quantum 
field is much less than the number-density of atoms, (a cc ) 
is always much smaller than unity. Therefore the mean 
intensity of the quantum field remains small compared 
to that of the control field even when the latter is turned 
to zero. 

C. decay of Raman coherence 

In the ideal scenario considered in sec. IV we disre- 
garded all processes resulting in a decay of coherences 
between the metastable states b and c. If such a de- 
cay at a rate 70 is taken into account, the group veloc- 
ity in the EIT medium cannot reach zero. Nevertheless 
an (almost) complete transfer of the quantum state of a 
propagating light pulse to a stationary matter excitation 
is possible. Introducing phenomenogical decay into the 
equation of the Raman coherence a bc results in a modi- 
fication of equation (|24|) according to 

?faa = ~-^ clAkz (jf t °bc + lQO b )j + F bc . (84) 

For formal consistency a Langevin noise term associated 
with the decoherence needed to be introduced as well. 
Such a dephasing trivially modifies the propagation law 
of the polariton 

i*(z,t) = exp^-70 J dr sin 2 9(rfj x 



x V^z-cjdT cos 2 9(t), oj +Fq,. (85) 

As expected dephasing simply results in a decay of the 
polartion amplitude (9) but does not prevent light stop- 
ping. A more detailed discussion of various realistic deco- 
herence mechanisms and their influence on the fidelity of 
the quantum memory will be presented in a subsequent 
publication. 

D. atomic motion 

We have shown above that adiabatically varying the 
envelope of the control field in time a light pulse of 
the probe field can be brought to a full stop. A spatially 
homogeneous variation of the group velocity can be real- 
ized, for instance, in the case when the control and probe 
field are propagating in orthogonal directions. However, 
in such a case the spin-wave component of the polariton 
has a phase, which is rapidly oscillating in z, since then 
Afc = kd — k\ = kd(l — cos 9) — kd- For an atom at 
position Zj one would have 

tfd*) = ItfMl expji^Zj jexpj-i^^} , (86) 

where a transformation from the rotating frame back 
into the lab frame was applied. While the second term 
corresponds to spatial oscillations with the small beat- 
frequency between pump and drive, the first term oscil- 
lates with an optical frequency. As a consequence the 
the polariton state is highly sensitive to variations in 
the atomic positions, or in other words would dephase 
rapidly due to atomic motion. To retain a high fidelity 
of the quantum memory it would be necessary to con- 
fine the motion of the atoms during the storage time to 
within a fraction of an optical wavelength, which is a very 
stringent condition. 

To avoid this problem nearly co-propagating beams 
can be used, such that Afc <C kd- In this case almost no 
photonic momentum is transferred to the atoms and the 
two-photon transition frequency would experience only a 
very small Doppler-shift. Such a configuration is there- 
fore robust with respect to atomic motion, as electronic 
and motional degrees of freedom are completely decou- 
pled. 

In a Doppler-free configuration propagation effects of 
the control field need to be considered. In the case when 
the probe field is weak compared to the control field at all 
times, the latter propagates almost as in free space and 
Q(z,t) — fl(t — z/c). Due to the resulting z-dependence 
of the group velocity, the trapping mechanism does not 
exactly preserve the shape of the photon wavepacket. 
However, when the probe pulse enters the medium with 
a group velocity much smaller than the velocity of the 
control field v° g <C c, the retardation of the control-field 
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amplitude can be ignored and fi(t — z/c) « Q(t). This is 
a good approximation as long as the time variation of the 
control field is sufficiently slow such that its spatial vari- 
ation across the probe pulse length is small. This implies 
that the characteristic time of change of the control field 
should obey T 3> L p /c, which is again easily satisfied. 

Before concluding we note that even in the case when 
no photon momentum is transfcred to atoms and there- 
fore atomic motion does not result in altering the phase 
of the spin coherence, in practice it is still desirable to 
suppress the motion. This is because the light beams 
normally have a finite cross-section and the atomic coher- 
ence is localized in the longitudinal direction and motion 
would tend to spread the localized excitation over the 
entire volume occupied by the gas. In this case the infor- 
mation about the pulse shape and the mode function will 
be lost after a sufficiently long storage interval. That is 
the reason why techniques involving cold trapped atoms 
and buffer gas were used to effectively slow the atomic 
drifts in recent experiments |2fj| , pT[| . 

VI. SUMMARY 

In conclusion we introduced the basic idea of quantum 
memory for light based on dark-state polaritons and dis- 
cussed their properties. We have considered the influ- 
ence of the main realistic imperfections on trapping and 
re-acceleration and have shown that the technique is ex- 
tremely robust. In particular, we demonstrated that an 
essential mechanism that enables the efficient quantum 
memory operation - the adiabatic following in polaritons - 
can take place despite of transparency window narrowing. 
Subsequent papers will discuss the influence of other de- 
coherence mechanisms on photon trapping and retrieval 
and will present the results of detail numerical simula- 
tions of the trapping procedure. 
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